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APPLICATION OF TWO-DIMENSIONAL VORTEX THEORY TO THE 
PREDICTION OF FLOW FIELDS BEHIND WINGS OF WING-BODY 
COMBINATIONS AT SUBSONIC AND SUPERSONIC SPEEDS 


By Arthur Wm. Rogers 
SUMMARY 


A theoretical investigation has been made of a general method for 
predicting the flow field behind the wings of plane and cruciform wing 
and body combinations at transonic or supersonic speeds and slender con- 
figurations at subsonic speeds. The wing trailing-vortex wake is repre- 
sented initially by line vortices distributed to approximate the spanvise 
distribution of circulation along the trailing edge of the exposed wing 
panels. The afterbody is represented by corresponding image vortices 
within the body. Two-dimensional line-vortex theory is then used to 
compute the induced velocities at each vortex and the resulting displace- 
ment of each vortex is determined by means of a numerical stepwise inte- 
gration procedure. The method was applied to the calculation of the 
position of the vortex wake and the estimation of downwash at chosen tail 
locations behind triangular-wing and cylindrical-body combinations at 
supersonic speeds. The effects of such geometric parameters as aspect 
ratio, angle of attack and incidence, ratio of body radius to wing semi- 
span, and angle of bank on the vortex wake behind wings of wing-body 
combinations were studied. The relative importance of wing vortices, 
the corresponding image vortices within tbe body, &nd body crossflow in 
determining the total downwash was assessed at a possible tail location. 


It was found that the line-vortex method of this report permitted 
the calculation of vortex paths behind wings of wing-body combinations 
with reasonable facility and accuracy. Α calculated sample wake shape 
agreed qualitatively with one observed experimentally, and sample results 
of the line-vortex method compared well with an available exact crossflow- 
plane solution. An empirical formula was derived to estimate the number 
of vortices required per wing panel for a satisfactory computation of 
downvash at tail locations. It was found that the shape of the vortex 
wake and the ultimate number of rolled-up vortices behind a wing depend 
on the circulation distribution along the wing trailing edge. For the 
low-aspect-ratio plane wing and body combinations considered, it appeared 
that downwash at horizontal tail locations is largely determined except 
near the tail-body juncture by the wing vortices alone for small ratios 
of body radius to wing semispan, and by the ος ii alone for large 
veer of that ratio. 
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INTRODUCTION 


Satisfactory aerodynamic design of high-Speed aircraft requires 
knowledge of the interference flow field resulting from wing-body-tail 
interaction. The behavior of the wing vortex wake in the presence of ' 
the body directly affects the air stream flowing past the tail surfaces. 
In particular, for certain relative sizes &nd positions of wing, body, 
and tail, the wing-body vortex wake produces large stream angles at the 
tail surfaces and nonlinear variation with angle of attack of the pitch- 
ing moment contributed by the tail. Therefore, it is desirable to be 
able to predict flow fields at the tail location for a given configura- 
tion in order to evaluate stability and control requirements. 


Furthermore, vapor-screen studies in supersonic wind tunnels have 
shown that the configuration of the vortex wake from the wings of some 
wing-body combinations differs markedly from the usual conception of & 
sheet whose side edges curl over to form a single pair of vortices. 
For example, figure l(a) shows the S-shaped cross section of a vortex 
sheet from each panel of a high-aspect-ratio supersonic wing and body 
combination at angle of attack.+ Here the ultimate vortex pattern far 
downstream of the wing consists of two vortices from each wing panel, 
both rotating in the same direction. Such & vortex wake results in & 
flow field different from the more well-known pattern, and the load on 
a tail situated in that flow field differs from the tail load associated 
with a single-vortex wake. It is important, therefore, to determine 
the conditions under which this unusual wake pattern occurs. 


Considerable investigation (refs. 2-12) has been devoted to the pre- 
diction of flow flelds behind wings, both plane and cruciform. Refer- 
ence 2 uses linearized conical-flow theory to calculate sidewash and 
downwash for a flat vortex sheet at the wing trailing edge and in the 
Trefftz plane, references 3 and 4 use supersonic potential doublet dis- 
tributions to determine downwash in the plane of a Plat wake and in the 
vertical plane of symmetry for any distance downstream of a wing. Ref- 
erence 5 employs pressure doublets to obtain general expressions for 
the induced velocities in space behind a wing, and references 6 through 8 
utilize line-yortex theory to predict sidewash and downwash in the flat- 
vortex wake and in space. The flow in the Trefftz plane behind super- 
sonic wings is treated in reference 9 by means of linearized conical-flow 
theory. The downwash based on a flat deflected vortex sheet and the 
downwash based on two deflected vortex lines are calculated in refer- 
ence 10. The rolling up of the trailing-vortex sheet behind wings is 
analyzed in reference 11, and reference 12 considers the motion of single 
vortices from each panel of a cruciform wing. The literature on the 
vortex wake behind wings of wing-body combinations, however, is still 
relatively meager (refs. 13 and 14). For wing-body combinations, ref- 
erence 13 accounts for the effect of the fuselage on the flow field by 


1A description of the vapor-screen technique is given in reference l. 
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considering the motion of two fully rolled-up vortices in the presence 
of a circular cylinder. Reference 14 treats a specific configuration by 
slender-body theory with the two assumptions of either a flat vortex 
sheet or two fully rolled-up vortices. 


The general problem to be considered here is the determinetion of 
the interference flow field behind wings of wing-body combinations 
(e.g., fig. 2) at subsonic or supersonic speeds. The dual purposes of 
this report are to present in detail and evaluate the method outlined 
in reference 13, and to apply it to the calculation of vortex paths and 
to the effects of wing aspect ratio, span loading, angle of attack, ratio 
of body diameter to wing span, angle of bank, and wing incidence on vor- 
tex paths and downwash at a possible tail location. 


Part I of this report presents the results of the application of 
the method. All calculations of span loading were made for triangular- 
wing and cylindrical-body combinations at a Mach number of 2.0. Wing 
aspect ratios of 2/3, 2, and ! are considered, combined with bodies 
whose radii are 0.2, 0.4, 0.6, and 0.72 times the wing semispan. The 
tail location selected for downwash computations in this part represents 
a missile-type configuration, for which the horizontal tail is assumed 
to be located in the body diametral plane 10 body radii downstream of 
the wing trailing edge. 


Part II of this report contains a detailed description of the line- 
vortex method. Απ examination is made of the effect on downwash at 
chosen tail locations of the number of vortices used. to represent the 
trailing vorticity. The nature of the stepwise integration method used 
is discussed, &nd solutions obtained by the stepwise integration method 
are compared with known exact solutions. Downwash is also computed at 
the tail of a high-tail airplane-type configuration, characterized by a 
horizontal tail 3 body radii downstream of the wing trailing edge and 
2 body radii above the body axis. 


SYMBOLS 
A aspect ratio 
c local wing chord 
Cr wing root chord in plane of symmetry 
ds directed line segment of a contour 
E elliptic integral of second kind 


Μεν unit vectors in x,y,z directions 


—— —— — — — — 
--------- y m A se _— 


Ap 


(| 


Uu,V,W 


r(y) 
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aspect-ratio correction factor defined in Appendix B 
lift 

Mach number 

number of vortices in summation 

Static pressure 


pressure differential across wing surface, P, - Pu 


PV? 


free-stream dynamic pressure, 


velocity vector, ui + vj + wk 

body radius 

local horizontal-wing semispan 

maximum horizontal-wing- semispan 

Spanwise position of wing vortex 

local vertical-wing semispan 

maximum vertical-wing semispan 

velocities in x,y,z directions (See fig. 2.) 
free-stream velocity 

chordwise distance from leading edge 


Cartesian coordinates of wind axes (See fig. 2.) 


load coefficient 


angle of attack of body axis, radians 

JW - 1 

angle of yaw, a sin Φ, radians 

spanwise distribution of circulation (defined by eq. (A1)) 


increment 
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€ downwash angle, E 
g — tan”? z 
A wing leading-edge sweep angle, 909 - w, radians 
HAN dimensionless rectangular coordinates of wind axes, Zi, 
(See fig. 2.) 
p air density of free stream 
σ,τ coordinates of source-point in xy plane 
p angle of bank, radians 
9179579, 
veloeity potentials 
Pos e. 
ω — S angle, tan E, radians 
Subscripts 
B body 
1,3 indices of summation 
D lower 
n maximum 
u upper 
W - wing 
LT linear theory 


SBT slender-body theory 
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l.- APPLICATION OF METHOD TO THE PREDICTION OF VORTEX PATHS 
AND DOWNWASH BEHIND TRIANGULAR-WING AND CYLINDRICAL-BODY 
COMBINATTONS AT SUPERSONIC SPEED 


General Description of Method 


The method used herein follows that of reference 13. The circula- 
tion distribution along the wing trailing edge of a given wing-body con- 
figuration is required as the initial condition. In this report, the 
circulation distribution is determined by modified slender-body theory. 
The circulation distribution is then replaced by & finite number of vor- 
tices which trail downstream and represent the wake vorticity. The 
effect of the afterbody is &ccounted for by vortices placed within the 
body &t the image position of each of the wake vortices. It is assumed, 
in general, that the flow changes in the stream direction are such that 
in the velocity potential equation 


the first term is negligibly small, that is, 


(1 - M?) Pay = 0 


so that 


. Thus the solution is independent of Mach number. Such an assumption 
of two-dimensionality is valid for slender, pointed wings and bodies at 
subsonic speed, and at supersonic speed when the entire plan form lies 
within the body nose Mach cone, or near M = 1 for more general con- 
figurations.? 


Consistent with this slender-body theory &ssumption, tbe bound vor- 
tices within the wing are neglected in computing vortex paths and down- 
wash. The trailing vortices are extended to infinity upstream and down- 
stream, and the induced velocities in crossflow planes are calculated by 
two-dimensional line-vortex theory and body potential crossflow. In the 
crossflow plane, therefore, the analysis considers the following picture: 


“See reference 15 for a discussion of these limitations. 
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| - Wing vortex repre- 
senting circula- 
2 Ε΄ tion distribution 


"EPA 
ARI à 
Body Image vortex 


un c Potential 


crossflow 


Sketch (a) 


The motion of each of the vortices in the crossflow plane is computed, 
end the results applied to the wing-body problem by relating time in the 
crossflow plane motion to distance downstream of the wing trailing edge. 


The wing and image trailing vortex system used in the stepwise cal- 
culations constitutes the only vortex pattern considered here. It is 
known that at high angles of attack and for large ratios of body diameter 
to wing span, viscous crossflow produces vortices above the body (e.g., 
fig. l and ref. 1) which significantly affect the flow field. A thorough 
understanding of the mechanism of formation, strength, and stability of 
crossflow vortices is &s yet lacking. Consequently, no attempt has been 
made to account for them. 


On the basis of & large number of computations made during this 
investigation, the following analysis was made of the effect of geometric 
parameters on the vortex wake and the induced flow field behind wings of 
wing-body combinations. Although only triangular wing and body combina- 
tions were studied, the general conclusions should be applicable to any 
wing plan form for which the circulation distribution is similar to the 
types contained herein. 
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Effect of Reduced Aspect Ratio, BA, and Span Loading 
on the Rolling Up of the Vortex Wake 


Figure 3 shows the calculated detailed rolling up of 10 vortices 
representing a typical wake behind a low-aspect-ratio configuration. 
Vortex coordinates for this figure are listed in table I. It can be 
seen that at 10 body radii downstream of the wing trailing edge, “Y of 
the 10 vortices have already spiraled into a single group which trails 
downstream in approximately the free-stream direction. The circulation 
distribution along the trailing edge of the subsonic-leading-edge wing 
of this configuration is of almost elliptic shape. The vortex sheet 
behind a wing of a high-aspect-ratio combination, represented in fig- 
ure 4, distorts much slower toward its final configuration and shows a 
different pattern. Vortex coordinates for this figure are listed in 
table II. The different shape of the vortex sheet is due to the change 
in span loading with reduced aspect ratio ΒΑ. 


The difference in span loading for the high- and low-aspect-ratio 
configuration is a consequence of the supersonic or subsonic leading 
edge, respectively. As a consequence of its supersonic leading edge, 
the wing trailing edge of the high-aspect-ratio combination of figure } 
has a theoretical circulation distribution which is linear from about 
the mid-semispan to the tip. This linear distribution is represented 
by the eight equal, uniformly spaced vortices shown, the remaining two 
vortices arising from the slender-body type of loading (eq. (3) of 
part II) inboard of the intersection of leading-edge Mach line and the 
wing trailing edge. Now, it is characteristic of a uniform distribution 
of isolated vortices along & line, such as those shown in figure 5, that 
the rolling up proceeds in the form of a symmetrical, S-shaped curve 
rotating &bout the centroid. The final configuration of the vortices in 
figure > will consist of two equal vortex cores, each containing four 
vortices, rotating symmetrically about the fixed centroid. In figure |, 
therefore, the S-shaped rolling up occurs for the uniform portion of ihe 
sheet, although the symmetry about and fixity of the centroid in space 
does not occur because the flow field is due not only to the 8 vortices 
but also to the other 12 wing vortices, 20 image vortices, and the body 
crossflow. The inboard portion of the vortex sheet displaces downward 
and outward in the conventional manner of figure 3. In summary, then, 
the effect of aspect ratio and the consequent change in circulation dis- 
tribution is to change the shape of the rolling-up vortex sheet as in 
sketch (b). The final vortex pattern for the high-aspect-ratio case 
depends on the extent of the span over which the circulation distribu- 
tion along the wing trailing edge is linear. A comparison of the calcu- 
lated S-shaped sheet and an experimentally observed vortex sheet is 
shown in figures 1(a) and 1(b). 
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Low BA High BA. Very high BA 
Trailing-edge circulation distribution 


Shape of rolling-up vortex sheet 


e, ; 


Final vortex pattern 
. Sketch (b) 


Effects of Geometric Parameters on Vortex Paths 
and Downwash &t & Tail Location 


In the following sections, calculated vortex paths and downwash are 
presented for a wing-body-tail combination at a Mach number of 2.0. The 
tail location at which downwash is determined is 10 body radii downstream 
of the wing trailing edge, a possible missile-type configuration. The 
vorticity shed from each wing panel is represented by a single vortex. 
For the comparison of the relative importance on downwash at the tail of 
the wing vortex, image vortex, and body crossflow, a single vortex for 
each wing panel suffices, in general, although the computation of the 
actual magnitude of downwash requires several vortices. In part II of 
this report an investigation is made of the effect of the number of vor- 
tices used to represent the trailing vorticity on the computed downwash 
&t chosen tail locations. 


Effect of aspect ratio.- Figures 6(a), (b), and (c) present the 
effect of aspect ratio on downwash, calculated at a missile-type tail 
location, for a small body and wing combination (r/sm = 0.2) at an angle 
of attack of 59. For the low-aspect-r&tio combination, the approximation 
that the total downwash is caused by wing vortices alone is good except 
in close proximity to the tail-body juncture. For increasing aspect 
ratios the approximation is good only &t greater spanwise distances from 
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the body. The reason for this result is simply that the downwash from 
the wing vortices decreases with increasing aspect ratio for the condi- 
tion of a given ratio of body radius to wing semispan. The wing lift 
decreases with increasing aspect ratio because of span load changes, and 
the effect of the body upwash therefore becomes more pronounced.  Slender- 
body theory predicts that the wake downwash angle exactly cancels the 
flow angle of attack inboard of the tip. It is interesting to note that 
inboard of the vortex location for the aspect ratio 2/3 wing (A = 4) the 
average downwash angle is roughly the negative of the angle of attack of 
the wing-body configuration. 


Effect of angle of attack.- Vortex paths behind wings of two 
triangular-wing and body combinations at various angles of attack are 
presented in figure T. It is seen that increased angle of attack results 
in & more pronounced inward and downward motion of the wing vortex with 
increasing distance u downstream, although for low-aspect-ratio com- 
binations the initial vertical motion is upward due to body upwash. 
Farther downstream the distance between the body and the vortex wake 
increases, with & corresponding decrease in the effect of body upwash. 
The vortices from the wing panels of the low-aspect-ratio configuration 
then move downward more rapidly because of their nearness to each other. 


Examination of the equations of motion of the vortex wake, presented 
in part II of this report, shows that for a given wing-body combination 
a single set of computations can be made Por.all angles of attack a. 

The results of figures 7(a) and (b) are replotted against ua in fig- 
ure 7(c), which then applies to all angles of attack, at the same Mach 
number. 


The effect of angle of attack on downwash is illustrated by fig- 
ures 8(a), (b), and (c) for an aspect ratio 2/3 wing and large-body com- 
bination, and by figures 8(d), (e), and (f) for an aspect ratio 4 wing 
and small-body combination. For the low-aspect-ratio configuration, it 
appears that the largest portion of the total downwash is contributed by 
the body upwash because of the large body. The difference in the shape 
of the total downwash curves of figures 8(a), (b), and (c) is due to the 
fact that the vortex wake passes farther above the horizontal tail as 
the angle of attack increases. For the high-aspect-ratio configuration 
the total downwash is not primarily caused by any single component 
because of the small body size. Wing and image vortices together with 
body upwash must be considered in calculating total downwash at any 
angle of attack for such configurations. 


Effect of ratio of body radius to wing semispan.- Another parameter 


governing the behavior of the vortex wake is the ratio of the body radius 
to the wing semispan, r/sm. Figure 9 presents the results for an aspect 
ratio 2/3 wing of fixed span, alone and in combination with bodies of 
different radii. The variation of the vortex strength with body size, 
given by equation (3) with y =r, is plotted in figure 10. The vortex 
strength goes to zero for r/sm = 1 since the entire wing then is 
enclosed by the body. It is seen in figure 9 that the initial slope of 
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the vertical displacement of the vortex path near the wing trailing edge 
increases in magnitude as r/sm increases. This is the result of 
increasing upwash in the body crossflow field coupled with correspond- 
ingly decreasing wing vortex strength (see fig. 10). The final uniform 
downward motion of the vortex pair begins when the body is far removed 
from the vortices, and the path lies above that for the wing alone case 
in which the vortex pair moves linearly downward from the trailing edge. 
The lateral motion of the trailing vortex is shown in the lower part of 
figure 9. Together with the decreased strength of the wing vortex for 
large values of r/sm, there is an initial outboard shift of the vortex 
at the trailing edge. For the isolated elliptically loaded wing, y/Sm 
equals x/! and it approaches 1 (the initial vortex location moves 
toward the wing tip) as r/sm increases from O to l. The lateral vor- 
tex motion downstream is affected by the sidewash component of the body 
crossflow field. The larger the body for & given wing, the more rapid 
is the inboard motion of the vortex. 


Figures 6(a), 11(a), 8(a), and 11(b) show the effect of r/sm on 
downwash for ratios of body radius to wing semispan of 0.2, 0.4, 0.6, 
and 0.72, respectively. The results shown in the figures can be antici- 
pated qualitatively. For an r/sm of 0.2, figure 6(a) shows that the 
downwash contributed by the image vortices almost entirely cancels the 
body upwash, except in close proximity to the body-tail juncture. Thus, 
for small ratios of body radius to wing semispan, downwash at the tail 
location is mostly given by the downwash caused by the wing vortices | 
alone except near the juncture. Figure ll(a), for an r/sm of 0.4, 
exhibits & decrease in downwasb from the wing vortices, and the total 
downwash is not as well given by the wing vortices as in figure 6(a) 
especially near the body-tail juncture. Figure 8(a), for an r/sm of 
0.6, shows that except near the body-tail juncture the downwash in this 
case is largely that given by the body upwash. Finally, figure 11(b) 
for an r/8m of 0.72, showing the same trend as figure (a) and fig- 
ure 8(a), indicates that for very large ratios of body radius to wing 
semispan, such as for canard configurations, the downwash at the tail is 
almost entirely determined by the body upwash. 


Comparison of Rolling Up of Vortex Sheet Behind 
Wing Alone and Wing-Body Combination 


Plane wing and body.- It is interesting to compare the rolling-up 
process of the vortex sheet behind a wing alone with that behind a wing- 
body combination under the conditions of sketch (c). 
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Elliptical Sm — 


Sketch (ο) . 


Both wings have the same aspect ratio, root chord, and Το. However, 
this maximum I occurs at the center line for the wing (with elliptical 
distribution) and at the wing-body juncture for the wing-body combination 
(with near-elliptical circulation distribution). Figure > of refer- 
ence. ll presents a picture of the rolling up of a vortex sheet behind an 
elliptically loaded wing as calculated by Westwater (ref. 16) in terms 
of the ratio of downstream distance d to rolling-up distance e. This 
ratio is related to H in figure 3 as follows (using eq. (5) of ref. 11): 


-($)e(8)(9) 
$- (2)@(4) 


d 
5 * 


alo 
| 
ala 


or 


where 
b wing span 
Cy, total lift coefficient 


K constant based on the shape of the trailing-edge circulation distri- 
bution 
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The value of K is 0.28 for elliptic loading. Now equation (10) of 
reference 11 is: 


C 
Io = sv, = 


where Ig corresponds to elliptic loading for the wing alone. When the 
last two equations are combined, 


Es (2) (0.28) RAS 


or 
a _ 2.2% i) S 
ER π \e/To/Vo 
In the present notation, 
a. 
S Em 


The lift coefficient of a triangular-wing and infinite cylindrical-body 
combination in terms of the wing alone lift coefficient is given by 
equation (38) of reference 15 as: 


p2 Y 
(Τη. = Cty ( 7 =) 


It is clear that while the same I, .is chosen in this comparison, 
the total lift of the combination is less than that of the wing alone 
because of the loss of exposed wing area. For the calculations of fig- 
ure 3, Doo equals 0.575 and Sm equals 3.75. These values being 
assumed in the equation for d/s, 


& - (338) - «e (3) 


For figure 3: 


Therefore the streamwise stations in figure 3 correspond to 0.2, 
O.4, 0.6, . . . 2 wing semispans downstream. The values of d/e in 
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figure > of reference 11 correspond to values of x/sm of 0, 0.233, 
0.977, 1.953, 4.370, and 8.040. Cross sections through the vortex sheet 
at these six stations have been drawn in figure 12 as well as the nearest 
corresponding sheet configuration from figure 3. From figure 12 it can 
be seen that the rolling-up process of the vortex sheet behind a 
triangular-wing and small-body combination closely resembles that for a 
wing alone under the conditions of the preceding sketch. 


Banked wing and body.- The qualitative difference between the 
rolling-up motion of a vortex wake behind an isolated wing and wing-body 
combination both banked at a given small angle and at angle of attack 
can be reasoned simply with the aid of sketch (4), illustrating condi- 
tions immediately behind the trailing edge. 


£ 

« 
κ 

/ 


SA 


Sketch (d) 


For the isolated wing, each panel vortex induces a velocity V, on 
the opposite vortex and both moye in the dotted direction with uniform 
velocity. The motion in this case is symmetric with respect to the n'u 
plane. Now the addition of a body, with potential crossflow symmetrical 
respect to the 7H plane, adds the same induced velocity V, at each 
vortex. The wing vortices in this case then move initially in the direc- 
tion of the resultant velocity Vp. Henceforth, the velocity due to the 
body crossflow is different at each vortex. Thus it can be seen that 
the wake motion is different from that behind an isolated wing, that is, 
completely nonsymmetrical. Similar reasoning for wing vortices of 
unequal.strength leads to the same conclusion. 


For a cruciform wing and body combination at 45° angle of bank and 
at angle of attack, the superposed velocity fields (from the wing and 
body) are both symmetrical with respect to the np plane. Therefore 
the computed vortex wake motion should agree qualitatively with results 
for isolated cruciform wings. 


-ma --ᾱν--- ----- ------ 
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II.- PRESENTATION AND EVALUATION OF METHOD 


The remainder of this report will examine in detail the line-vortex 
method outlined in reference 13. The manner by which the initially 
required circulation distribution was obtained is first discussed. The 
method of replacing the wake circulation distribution and the afterbody 
by discrete vortices then follows. The determination of the flow field 
in transverse planes and the stepwise integration technique are next 
explained. After a few sample calculations, an evaluation is made of 
the errors inherent in the stepwise integration and of the effect of the 
number of vortices used on the computed downwash distribution at chosen 
tail locations. 


Circulation Distribution at Wing Trailing Edge 


Before the stepwise calculation of vortex paths behind the wings of 
wing-body combinations can be started, it is necessary to know the cir- 
culation distribution along the wing trailing edge. In general, the 
circulation distribution T(y) and the wing span loading are not equiva- 
lent. The conditions for equivalence are derived in Appendix A. The 
method of this report has been applied to triangular-wing configurations 
which have readily obtainable circulation distributions. However, the 
stepwise calculation method is applicable to configurations with any 
wing plan form. 


For a plane-wing and cylindrical-body combination, the span loading 
15 equivalent to the circulation distribution and can be obtained by 
chordwise integration of the complete pressure coefficient (Ap/a)y- 
Equations for (Ap/a),, are given in the slender-body analysis of refer- 
ence 15. For such configurations, the circulation about the wing at any 
spanwise station y is given by linear theory as: 


I(y) "A dx; (1) 


The notation for this and the following equations is defined in fig- 
ure 13, which shows the more general case of & noncylindrical body. 


Equation (11) of reference 15 furnishes the wing loading coefficient 
for a plane-wing and cylindrical-body combination as: 


as (1-2) 
(2 mu με BE 
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Substitution of equation (2) in equation (1) and integration yields: 


=o V (n? sg) e y se (3) 


When the body is not cylindrical, Bquetion (2) is modified to: 


(2) E HERE as - =; 


The last term of tión (4) is due to the spanwise velocity and 
does not contribute to the circulation distribution. Hence one obtains 
the required circulation distribution by insertion of only the first 
term of equation (4) in equation (1) and integration. The result is 
identically equation (3). This shows that the circulation distribution 
I(y) along the trailing edge of a plane wing on a noncylindrical body is 
dependent only on the cross section of the configuration at the trailing 
edge, as could be anticipated. The circulation is, in fact, the jump in 
the crossflow potential at the trailing edge and hence must agree with 
the slender-body result, equation (3). 


For cruciform wings on a cylindrical body at angle of attack and 
sideslip, the complete velocity potential (eq. (47) of ref. 15) is: 


® = sleet]. (a4 r^ |r.? cos ο. z) l+ ος E. 
E p ο; τα - 3 
4 2 Ἐπ j τι 
ΑΝΘ Lol. 
4 4 2 
lr^ cos? 90 + s* (1 + 27) - 855 c =) (2+ n cos 5ο } x 
S τι sf 
V. p! 4 rt 
Voatz + — πο. r? cos ο EF) | 
y E i 4 
2 PI t 
2 


4 N2 i 
tn l- => + duró cos? 20 + tt (15 23 + 
E : τ 
4 4 >94 
r r 2 |Z 
ete € + E € E) I COS 20 P) + VoP'y (5) 
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where 


r4? = y? + z?, O = tan”? 


«dits 


and 


+ means O«0-«mx 


- means x <@ < 2x 


The transformation to a configuration pitched a radians and banked 
Q radians is for small angles a: 


at =a cos Y, Bt = a sino (6) 


Now for banked cruciform configurations, the span load distribution 
is not equivalent to the circulation distribution. Hence, instead of 
equation (1) one must use the general relation between TI and the jump 
in the velocity potential at the wing trailing edge: 


i (7) 


For horizontal and vertical surfaces respectively, one obtains by 
substitution of equation (5) in equation (7): 


ry) «το NA - τ (8) 


Sm 


T a (tm^ - z?)(tm^z? - r4) (9) 
m2 


r(z) 


It is thus seen, by comparing equations (8) and (9) with (3), that 
the circulation distribution for a cruciform wing and body combination 
is derivable by assuming that each wing acts independently of the other 
as part of a plane wing and body combination at an angle of attack given 
by equation (6). This can be seen also by the linear superposition of 
potentials in equation (5). 


The foregoing equations are based on the slender-body theory of 
reference 15, which is postulated for slender wing and body configura- 
tions at subsonic, transonic, and supersonie speeds. In order for the 
equations to apply at supersonic speeds, the entire configuration must 
lie well within the body nose Mach cone. Experimental data on models 
conforming to these restrictions agree well with the theory. In order 
to apply the results of slender-body theory to nonslender wing and body 
combinations, an aspect-ratio correction factor is employed. This factor, 
which for triangular-wing configurations is the ratio of the lift of a 
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triangular wing calculated by linear theory to the 11Pt calculated by 
Blender-body theory, is discussed fully in Appendix B. The results may 
be summarized as follows: For subsonic-leading-edge wings (B tan w <1), 
equations (3), (8), and (9) should be multiplied by the factor 


xl (1 Paro) 


For supersonic-leading-edge wings, (B tan w >1), these equations should 
be multiplied by equations (B7) or (B18) for the ranges of y indicated: 


k = TIU a (BT) 

ες ΠῚ m y 
for 

qoc ugue. 

B tan ω 
2 2 c 
k s fon +2) (em - eni A M 
(sm - y”) (B? tan? w - 1) B ten w(Sm - y) 


yB? tan? ω +Sm 
(sm + y) sin“? e | - (B18) 
Btan w(sm + y) 
for 
-r 
ες pq E 
umb B tan ω 


Replacement of Wake Circulation Distribution 
and Body by Discrete Vortices 


After the circulation distribution TI(y) along the wing trailing 
edge has been determined, the next step is the replacement of the conse- 
quent wake circulation function by a finite number of vortices. A plot 
of the circulation function, equation-(3), for example, looks like that 
portion of the solid curve above the wing in sketch (e): 
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Sketch (e) 


The loading over the body, which can be calculated from the corre- 
Sponding body pressure coefficient (eq. (12) of ref. 15), is not explic- 
itiy required here. It will be seen that this loading is automatically 
accounted for by the vortex model set up for the wing loading. The 
curve I(y) over the wing can be replaced by a suitable number of step 
functions (dashed lines in sketch), each of which corresponds to a vor- 
tex of strength equal to the step height and located at the step abscissa 
as shown. Although an infinite number of step functions would be 
required to duplicate the given curve, in practice the number selected 
are the fewest which render possible a reasonably accurate prediction of 
the trailing-vortéx sheet behavior. Single-, three-, and ten-step func- 
tion distributions have been used herein. Quantitative results will be 
presented later concerning the number of vortices used to replace the 
circulation distribution. For a single-step (one vortex) approximation, 
integrating equation (3) for the area under the TI(y) curve and equating 
the result to a single rectangular area lead to the simple result shown 
in figure ll, namely, that the spanwise distance from the body of a 
single vortex representing the circulation distribution on one panel of 
a plane, subsonic-leading-edge wing and body combination is approximately 
constant and equal to about 0.76 times the exposed wing semispan for all 
ratios of body radius to wing semispan. In this approximation by a 
single vortex the value of Imay at the wing-body juncture is assigned 
to the vortex. In all cases the two panels of the wing are treated inde- 
pendently. A circulation distribution which, unlike the preceding sketch, 
is asymmetrical across the wing trailing edge will, therefore, require 
an asymmetrical vortex distribution. 


The effect of & circular body in the presence of this wing vortex 
distribution is accounted for by placing an image vortex for each wing 
vortex within the: body at the inverse point, as indicated by sketch (f). 
The image vortex cancels the velocity normal to the body due to the wing 
vortex. 
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Sketch (f) 


The foregoing discussion applies to the trailing-vortex filaments 
from the wing and the image vortices contained within the body. These 
filaments can be considered joined by bound vortices within the wing 
and adjacent body to form horseshoe vortices as shown in sketch (g). 


Sketch (g) 


Consistent with the assumed two-dimensional nature of the solution, 
the segments of the vortices bound within the wing and adjacent body are 
neglected. It is then assumed that the remaining filaments extend 
upstream to infinity. 


The portion of the total lift of 

the wing-body combination carried 
by the wing and by the body is rep- 
resented by the bound filaments of 
length lo within the wing and 2, 
within the body, respectively. A 
plot of T(y) across the body is, 
therefore, obviated since the single 
vortex representing the body loading 
automatically appears at the inverse 
(image) point of the wing vortex. 
This vortex model, due to Lennertz 
(ref. 17), assumes that the ratio 
of body to wing lift is proportional 
to the ratio 1,/lo2 (11 + lo being 
the "effective semispan"), and is 
known from experiments to represent 
the distribution of lift between 
wing and body with good accuracy, 
deviating but slightly from slender- 
body theory. 


Determination of the Flow Field in Transverse Planes 


The nonlinear differential equations governing the downstream motion 
of the vortex wake require a stepwise solution for the vortex positions 
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at successive streamwise stations (or time increments). Consider a 
vortex pattern &t &ny downstream station such as section A-A in figure 2. 
(The image vortices inside the body are not shown in the figure.) It is 
assumed that every vortex influences the flow field according to the 
Biot-Savart velocity law of planar, incompressible-flow vortex theory, 

as illustrated in the sketch (h). This assumption will be justified 


later. 
+W " NV 
d abs 
AN 
τν ve os 
(Acn) 
"M -— 
λ 1 : 
e η} ! Sketch (h) 
The velocity V4 induced at the ith vortex by the jth vortex is: 
ES 
Vi = 10 
i" eub (10) 
The vertical component, wi, of Vy 15: 
r -À 
Np rns V EE Puma Y (11) 
exr (Az - Ay)” + (Mi - 4) 
The horizontal component, vy, of V4 is: 
° -I 74 - T 


ear (Ay - A? + (πὶ - 037? 


In general, for & flow field containing n vortices, the velocity 
components induced &t the position of the ith wing vortex by the other 
n-l vortices are: 


n T4(Aq - Ay) 
l J J 
ua) Og A4)? + (ng - 14)? = 
τ 17Aj 17 
1 
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Qo T.(n, - 93) 
ud Ro ua e 
Ji 


Equation (10) is derived from consideration of an infinitely long, 
incompressible-flow vortex filament, and is, therefore, applicable only 
to two-dimensional flow. For supersonic speeds, the use of equation (10) 
is consistent with the slender-body theory underlying this study for 
these reasons: Reference shows that an infinitely long vortex filament 
parallel to a supersonic stream obeys the Biot-Savart law of incompres- 
sible flow. Omission of the bound vortex and wing chord loading causes 
the difference between initial and asymptotic downwash at the wing trail- 
ing edge seen in figure 15 reproduced from reference 7; the downwash 
along the wake center line approaches the asymptotic value within a chord 
length behind the trailing edge for the triangular wings of PA< 1. 
Hence the asymptotic downwash given by a two-dimensional trailing-vortex 
system has been used. 


superimposed on the velocity field due to the wing and image vortex 
system is the body potential crossflow velocity field. From the poten- 
tial: 


Pi = Vg m (15) 


where z' is measured from the body axis (instead of the wind axis), one 
readily obtains the velocity components induced by the body: 


ὃ "—-—- 
wW = A = voar? A (16) 
dz (y? + z! ) 
OQ, . -2V,ar?yz! - 
diri pun ος (1T) 
dy — (y? + 21%) 
or in dimensionless form: 
2 2 
"NE d M (18) 
Yo (A? + 11%) 
- -2 t 
ds E (19) 
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where 
η! =97+|41tna=0"+ pa (20) 


This transformation from the usual potential equations is neces- 
sitated here by the fact that the coordinates are not body centered, and 
the body is inclined, with respect to the free-stream direction. 


These last equations assume that the body angle of &ttack is small 
enough so that tan a za, cos a € 1 (which is true for a up to approxi- 
mately 20°), and that the contribution of the crossflow to the local 
streamwise velocity is negligible. As stated previously, viscous cross- 
flow separation around the body, with its consequent vortex wake, is not 
considered here. 


Stepwise Determination of Vortex Paths 


The replacement of the wing trailing-edge circulation distribution 
by step functions and associated vortices, together with the placement 
of image vortices within the body, as detailed above, provide the start- 
ing point for the stepwise calculations to be described. The spanwise 
and heightwise coordinates of all vortices at μΞ0 are known. The 
downstream incremental motion of these vortices will next be determined. 


Now at any point in a transverse (yz) plane, the lateral velocity, v, 
is given by the sum of equations (14) and (19). The vertical velocity, W, 
is similarly given by the sum of.equations (13) and (18). The streamwise 
velocity, u, perpendicular to the plane, is everywhere V,- Consider 
the point where a vortex line passes through the transverse plane. In 
a time interyal At the vortex filament moves with the fluid a vertical 


distance Az = w At. Since At = Axo, Az = (w/Vo)^x, then nondimen- 
slonally: 
An = = Ap (21) 
ο 


Similarly the dimensionless lateral movement is: 


Ad Am (22) 


The last two equations then furnish points in transverse planes a 
distance AH apart, the locus of which is the path of a vortex filament 
moving in the wing-body flow field. With each vortex, & "strength" Iso 
(rather than Ty) will be associated in order that equations (13) 
and (14) will ος yield the velocity ratios w/V, and v/V,. 
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In the plane u= 0, which contains the wing trailing edge, as long 
as the wings have no incidence with respect to the body, all the vortex 
filaments have coordinates η =0. Therefore, the total spanwise veloc- 
ity induced at the location of any given wing vortex can readily be seen 
to be zero, from equations (1%) and (19). In fact, the initial motion 
of such a vortex sheet is always one of pure vertical displacement. The 
total vertical velocity summed from equations: (13) and (18) is multiplied 
by a suitable increment Ap yielding an incremental distance An moved 
by the given vortex line, according to equation (21). Values of An are 
thus calculated for all the wing vortices, the entire group of which is 
so transported to the plane pn - Ap. The coordinates of each wing vortex 
in this new plane are found simply by adding the computed increments to 
the coordinates &t the previous station. Of course, symmetrical proper- 
ties are used wherever possible to obviate the calculation of the paths 
of vortices from each half of a wing. Corresponding to the new location 
of each wing vortex line at the station u= An, the image vortices 
within the body are repositioned according to the following formulas: 


À 
Maga === (23) 
7 a (η + u tan a)? 


η + p tana 
"image = -H tan a + (24) 


A? + (η + p tan o)? 


where (λ,η) are the coordinates of the wing vortex at the new station 
H = Ap. This readjustment of the image vortex positions, required by 
the downstream displacement of the inclined body away from the free- 
Stream direction, is made because the image vortices are bound within 
the body; the image vortex displacement, in fact, measures the amount 
of lift carried by the afterbody (see ref. 13). 


With the new positions of wing and image vortices this procedure is 
repe&ted. That is, the total lateral and vertical velocity induced at 
each wing vortex location is calculated and multiplied by an incremental 
distance AH. The new wing vortex positions are obtained by adding the 
computed increments to the original coordinates. The image vortices are 
then repositioned according to the new wing vortex locations, and the 
stepwise calculation continued to, the desired downstream station. 


Ihis stepwise procedure is simply & method for integrating numeri- 
cally n simultaneous differential equations of motion of n/2 wing 
vortices for the downstream paths. A closed solution for the three- 
dimensional paths is, in general, extremely complicated, if not impos- 
sible. While the equations of the vortex motions in a transverse or 
crossflow plane can be written, introduction of time dependence (equiva- 
lent to the streamwise coordinate) renders a closed solution extremely 
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difficult. Reference 13, using complex variable notation, cites explic- 
itly the entire path of a pair of vortices in the presence of a cylinder 
in the transverse plane. The equations were derived by Villat in refer- 
ence 18. Even for this simple case, where n = 4, recourse is necessary 
to numerical solution for the time or streamwise coordinate. This solu- 
tion will be analyzed in greater detail subsequently, as well as the 
manner of choosing the proper spacing AH of successive stations. 


Sample Calculations 


Three examples will now be presented to illustrate the calculative 
procedure and resulting vortex paths. They are: 


l. A plane-wing and cylindrical-body combination at angle of 
attack. 


2. A cruciform-ving and cylindrical-body configuration at angle 
of attack and bank. 


3. A plane-wing and cylindrical-body combination with the body at 
“angle of attack and the wing panels at differential incidence. 


The procedure has been systematized by prepared computation forms 
such as table III which is used when the circulation distribution on one 
wing panel is replaced by & single vortex. 


Example 1, plane triangular-wing and cylindrical-body combination 
at angle of attack.- The data for this example are: wing aspect 


ratio - 2/3, ratio of body radius to wing semispan = 0.6, angle of 
attack = » and Mach number = 2.0. 


The leading edge of the wing is subsonic (B tan w= 0.2912 <1), so 
the trailing-edge circulation distribution Τ(Υ) is calculated from equa- 
tion (3) and multiplied by the aspect-ratio correction factor k = 0.917 
(see Appendix B). A plot of p(y) vs. y like that of figure ll is 
obtained, with a maximum ordinate I/V, = 0.12796 in this example. 
First, this distribution is replaced by & single vortex per wing panel, 
with the location consequently given by figure 14 as s! /(Sm - τ) = 0.763, 
or A, = 1.509 (r/sm = 0.6). The vertical coordinate of the vortex at 
the wing trailing edge (u = 0) is zero in this case. The dimensionless 
spanwise coordinate of the image vortex is, therefore, As = LA = 0.663. 


Referring now to table III(a), the above values of rA » Ay» η 
No ; and Mp are filled into the proper boxes of columns (2), (3) and . By 


symmetry, the values of 4, A,» Nar end η, can be readily filled in, and 
an increment Ap is chosen (1 here). Columns (5) through (9), based on 
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equations (13) and (14), are then computed in order, the last two giving 
the vertical and horizontal velocity components induced at vortex l by 
the other three vortices (columns and (9), respectively). Rows 
through 63, based on equations (18) through (22), next add the potential 
crossflow velocity components at vortex l. Thus, rows and 6 give 
the new coordinates of vortex 1 for use in table III(b), while the 
bottom calculation furnishes the coordinates of the corresponding new 
image vortex 2 position. 


Columns (8) through (1) of the table III(b), for station p = 1, can 
now be filled in and the procedure repeated, as illustrated. This step- 
wise computation is continued for as many stations downstream as desired, 
each part of this example representing positions at increments of one 
body radius (Au = 1) downstream of the wing trailing edge. The results 
are summarized graphically by figure 16. 


For the same given data, the calculated circulation distribution 
for one wing panel is now replaced by three vortices of equal strength 
EV = 0.0h2653. Their spanwise locations are determined graphically, 
equating areas under the curve of T(y) vs. y, and in this example are: 
Ay = 1.652, A, = 1.558, and λα = 1.316. The three image vortex loca- 
tions are again calculated as the inverse points, giving: A, = 0.75988, 
As = 0.64185, and Ag = 0.60533. The results of this case are shown in 


figure 17. 


le 2, cruciform triangular-wing and cylindrical-body configura- 
tion at angle of attack and bank.- The configuration data of the pre- 
ceding example are again assumed in addition to which another pair of 
wing panels now exists, forming a cruciform wing-body &rrangement banked 
45°, The stréngths of the assumed four vortices originating from the 
four wing panels are identical, as required by equations (8) ana (9). 
These equations with equation (6) provide the required circulation dis- 
tribution [(y) and [(z) which is corrected for aspect ratio as explained 
in Appendix B. Figure 1} again furnishes the initial positions of the 
four wing vortices, whose corresponding four image vortices are then 
located at the inverse points. The calculative procedure follows that 
described above. Because of planar (xz) symmetry, only the paths of a 
pair of vortices on one side of the body need be computed. The results 
of this case are presented in figure 18, which shows that there is a 
tendency toward "leapfrogging," sketch (1)9 although the process appears 
to be retarded initially by the body. 


SThe phenomenon of "leapfrogging" vortices is discussed in detail 
for isolated wings in reference 12. 
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"Leapfrogging" paths in crossflow plane 


Sketch (1) 


le 3, plane-wing and cylindric&l-bo combination with the 
body at angle of attack and the wing panels at differential incidence.- 
The given data' for this case are the same as for the first example with 
two exceptions. One wing panel is deflected and the other is at zero 
incidence with respect to the body axis such that the vortex from each 
panel has the same strength. Now the vortices from both wing panels 
are rotating in the same direction. The vortex shed Prom the deflected 
panel is assumed to start above the xy plane because of the angular 
displacement about the hinge line of that panel's trailing edge. It is 
further assumed that the vortex from each panel is shed at the same span- 
wise station. 


Figure 19 illustrates an effect characteristic of wings differen- 
tially inclined. The vortex from the wing panel which has no incidence 
to the body moves essentially the same way &s in the plane wing case. 
However, the vortex from the wing panel which is at a negative angle of 
attack travels sharply upward and inboard. This marked motion is due 
to two effects: First, the vortex starts at the trailing edge with some 
η >0 because of the negative angle of attack of the panel, and second, 
the body crossflow and image vortex-induced velocities are cumulative in 
this case instead of being subtractive as in the plane case. At & suf- 
ficiently high angle of attack, it may be possible for the rising vortex 
to jump over the body to the same side as the other vortex. Such an 
occurrence would cause a sharp increase in rolling moment from the tail 
surfaces. l 


For the case of a pair of wing panels inclined together, the trail- 
ing vortex paths can be estimated qualitatively by consideration of the 
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vortex picture in the crossflow plane. Thus, for a pair of wing panels 
inclined at a positiye angle of attack on a body whose axis is parallel 
to the free stream, the vortex paths trail below the free-stream direc- 
tion with concavity upwards since there is no body upwash and the effect 
of the image vortices is to depress the wing vortex wake. For a pair of 
wing panels at a negative angle of attack on a body at a positive angle 
of attack, the trailing vortices are inclined upward since the body 
crossflow and vortex-induced velocities are additive. 


These and all other computations except those with 10 vortices per 
wing panel were performed with a desk calculator. It was found that the 
solution of vortex paths downstream of the wing could be accomplished in 
about fifteen minutes per station for a single-wing vortex, and two hours 
per station for a three-wing vortex scheme. In general, the number of 
computations per solution increases approximately as the square of the 
number of vortices assigned per wing panel. 


Evaluation of the Effect of Stepwise Integration 
on Vortex Paths 


in the present line-vortex method, the principal factors governing 
the l&bor expended in the solution to & given problem are the size of 
the "time" increment Ap and the number of vortices used to replace the 
trailing vorticity. The first factor will now be discussed in detail. 


The necessity for a stepwise solution for vortex paths stems from 
the difficulty of integrating n simultaneous differential equations of 
vortex motion (eqs. (13) and (14), together with eqs. (18) and (19)). 

The differential equations are therefore solved with a small, finite Ap. 


In this report, the use of equations (21) and (22) implies the use of 
Euler's linear integration method, that is, 


7 avy 
aay = Ya, + ( dy ). T" 
where 


VQ = Ay Or My and i=l, 2, ds... ņ/2 


A Taylor series expansion shows the higher order terms neglected: 


um = Od, (πὸ | Ap + E), o rl = x sues 


Thus, it may be seen that the use of Euler's method without appre- 
ciable error requires either the second and higher derivatives to be 
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small or Au to be small. The following discussion indicates the impor- 
tance of the higher order terms for vortex motion. 


Comparison of known exact solution of vortex motion with stepwise 
integration.- A simple example of the mutual interaction of two vortices, 


for which the exact solution is known, furnishes an insight into the 
nature of the cumulative error incurred by this stepwise-integration 
method. It is easily shown by the methods of hydrodynamics (e.g., p. 320 
of ref. 19) that two vortices of strengths ΓΙ and Tp. will, due to their 
own mutual influence, each rotate about their common “centroid” in a cir- 
cular path with constant angular velocity which is equal to (Is + r,)/2x 
divided by the square of the distance between the vortices. This motion 
is illustrated in figure 20. The solid lines radiating from tbe centroid 
represent the ends of a constant time increment (proportional to Ap) and 
the solid, spiral-like curves are the corresponding paths computed by a 
stepwise approximation. The dashed radial lines and spirals relate to a 
time interval half of that used before. It appears from figure 20 that 
the error or discrepancy between circular and spiral paths increases 

while the rate of growth of this error decreases with time (or distance 
downstream of the initial position). Further, it is seen that a decreased 
size of time increment results in & decreased error, for this example, in 
approximately the same ratio. Of course, this case of completely circular 
vortex motion cannot be solved accurately by Euler's linear scheme without 
using extremely small increments. However, for many of the computed vor- 
tex paths behind the configurations considered in this report it will now 
be shown that Euler!s method is sufficiently accurate. 


Effect of size of increment Au on computed vortex paths.- The 
error incurred by using the linear integration method can be judged by 


calculating vortex paths with different increments Apu and extrapolating 
the results to AH = O. This has been done for single- and three-vortex 
schemes at a= 59 and typical results appear in figures 21 and 22. 
Both figures show that the stepwise method employing finite increments AH 
yields results which converge fairly rapidly to the exact solution (for 
which Δμ--»άμ--»0). Here the error, using a given value of AH, is 
indicated by the difference between the value of A or n extrapolated 
to Au-0 and the value of Aorm at the given value of Ap As in 
figure 20, the error is seen to increase with time or distance u down- 
Stream of the wing. For a single-vortex scheme (fig. 21), the error is 
approximately proportional to pAp and is small, in general. The error 
becomes greater and less predictable, for the same increment size, with 
a multiple-vortex scheme (fig. 22). A single-vortex scheme is therefore 
less subject to integration error than a multiple-vortex scheme. The 
reason why the error increases nonlinearly for the three-vortex scheme 
of figure 22 can be explained by reference to figure 22(c) which shows 
the projection of the vortex paths on the crossflow plane. For small 
values of μ the vortex paths are fairly linear, with the vortex posi- 
tion error given by figure 22(a) being small and essentially linear 

with AH. Further downstream, as the effect of the rolling-up process 
becomes marked, the paths are seen to become spiral-like. Each vortex- 
position error curve becomes incressingly nonlinear as the vortex path 
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approaches a maximum or minimum vertical position in the crossflow plane. 
The use of smaller increments Ap then reveals the tendency toward con- 
vergence of the path to one of a lesser radius of curvature, as illus- 
trated previously for the case of truly circular paths. 


From these results it appears that the selection of a suitable 
increment Ap for any given case should be made by trial calculation 
with several values of Ah for a few stations downstream of the wing. 
As long as the computed vortex paths are relatively linear in the cross- 
flow plane, the largest value of Ap should be chosen which permits the 
extrapolated error in vortex location to be the maximum tolerable ai the 
furthest downstream station at which downwash is desired; Ap should be 
decreased locally wherever the vortex p&ths appear to be approaching a 
maximum or minimum height in the crossflow plane. 


More accurate numerical integration schemes than Euler's method can 
be found in reference 20. Curve-fitting and extrapolation formulas are 
presented which by taking account of the higher order derivatives of the 
path enable one to maintain sufficient accuracy of vortex positions with 
a given Ap even when the paths are markedly nonlinear. 


Comparison of computed paths with exact solution for vortex paths 
in crossflow plane.- The exact solution known for the vortex path in the 


crossflow plane corresponding to the single-panel voriex representation 
mentioned previously can be used to check the accuracy of the stepwise 
integration. Equation (VI-h) of reference 13 cites the result obtained 
by Villat (ref. 18). In the present notation, 1t is: 


2a (1- &) AE: HA? (£2.-1)* «ος (ids ja E 2X5 (97-1) -1) 
ο ο 


lixVo (£2-1)? + 2 Nz el +1 
(25) 

where 

2 2 2 

ES = AT + (n + pa) 

λο m (^) =o 
At an infinite distance downstream the vortices are infinitely far 
removed from the body (η + pa—>e) so that the asymptotic spacing 
2 (A) use 9 the vortex pair is given by the relation: 
2 
r - 1 p θλο(λο” -1) 
2 (À - —— In 2(A = Dia LA Ya τ (26) 
Ίμ-ω 2xVo ) um o( xz) enV, Ao +1 
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Two solutions by the stepwise method have been checked against 
equations (25) and (26). This has been done by inserting pairs of values 
of the coordinates pu end η, from the stepwise-calculated paths, into 
equation (25) and solving for the A coordinate by trial. The compari- 


son appears in figure 23, showing the results to be practically coincident. 


Effect of Number of Vortices on Downwash at Tail Locations 


Comparisons were made of the downwash along a line, representing a 
horizontal tail plane, resulting from the separate presence of one, 
three, and ten vortices, all derived from a given circulation distribu- 
tion. Possible tail locations for both a missile-type and high-tail 
airplane-type configuration were chosen. 


. Missile-type configuration.- For the former, p = 10 was selected; 
that is, a tail location 10 body radii downstream of the pl a 
edge. Therefore downwash was first computed along the line = 10, 

n = - 10 tan a, 1<A< 6, which is a horizontal line ΟΕ gs body. 
center line of the missile-type configuration. 


Figures 24(a) and (b) show the effect of number of vortices on down- 
wash at the chosen missile tail location with and without components of 
the body potential crossflow in the stepwise computations. These omis- 
sions were made for two reasons: One was to determine the extent ‘to 
which the stepwise calculations could be simplified without obtaining 
dissimilar results in downwash, and the other was an attempt to account 
for the flow around the body at high angles of attack.. As an assumption 
for the high-angle-of-attack crossflow, the potential crossflow (eqs. (18) 
and (19)) was replaced by equation (18) alone with n' = 0. This assump- 
tion implies that the velocity distribution on the lee side of the body 
is the same at any vertical station as along the horizontal diameter, as 
indicated in sketch (j). 


P 


Potential crossflow. Assumed center line crossflow. 
Sketch (j) 


— -- - ----------- ----- 
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For the aspect ratio 2/3 configuration of figure 24(a), with 
r/sm = 0.6, the mmber of vortices used has a negligibly small effect on 
downwash except when the body crossflow is entirely omitted in the step- 
wise calculations. Observance of the various vortex positions at the 
tail location, shown at the top of figure 24(a), shows that the vortices 
from the wing panels are subject to a relatively strong body crossflow 
so that their computed positions are greatly in error if the body upwash 
is neglected. Thus, when at least the upwash component of the body cross- 
flow (the center. line crossflow in sketch (3)) is added to the stepwise 
calculations there is. no appreciable effect of either the number of vor- 
tices used or the strength distribution for a given multiple-vortex ` 
scheme. Figure 24(a) also shows that only a small change in downwash 
results from replacing the potential body crossflow by the assumed center 
line crossflow. | 


For the aspect ratio ! wing and body combination of figure 24(b), 
with r/sm = 0.2, it is seen that there is little effect of body cross- 
flow on the downwash at the tail. This could be expected because of 
the relatively large wing. While the use of & single vortex does not 
‘adequately approximate the magnitude of downwash over a tailspan equal 
to the wingspan, a 3-vortex scheme appears to be as satisfactory as a 
lO-vortex scheme. Reference to the computed vortex locations shown at 
the top of the figure indicates that the vortex sheet from a high-aspect- 
ratio wing at low angles of attack rolls up so‘slowly and extends over 
such & large spanwise distance that & single vortex cannot give the cor- 
rect distribution of downwash. 


At higher angles of attack, the vortex wake will be even further 
from the horizontal tail plane because the body will be inclined below 
the free-stream direction more than the vortex wake according to fig- 
ure T. Although the vortex strengths increase directly with angle of 
attack, the downwash at the tail decreases very rapidly with distance 
from the wake. Consequently, tbe effect of number of vortices should be 
smaller at higher angles of attack. 


Airplane-type configuration.- Now a possible horizontal tail loca- 
tion for a high-tail airplane-type configuration is u = 3 and 2 2, 
that is, 3 body radii downstream of the trailing edge of the wing and 
2 body radii above the body axis. Figure 25 presents the effect of the 
number of vortices used on the downwash at this tail location. The 
aspect ratio 2/3 configuration of figure 25(a), with r/sm = 0-6, shows 
no noticeable effect of number of vortices on either the distribution or 
magnitude of downwash simply because of the distance of the horizontal 
tail plane above the vortex wake. Neither is there any appreciable effect 
of assuming the body crossflow to be the center line crossflow. The 
aspect ratio 4 configuration of figure 25(b), with r/sm = 0.2, shows 
that 3 vortices give the same downwash as lO vortices. However, from 
figure T the vortex wake behind the wings of both high- and low-aspect- 
ratio configurations at small values of yu is seen to be fairly close 
to the free-stream direction at all angles of attack. Therefore it can 
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be expected that with increasing angle of attack, the tail plane 
approaches the wing vortex wake and the use of at least three vortices 
to represent the weke behind wings of such high-tail airplane-type con- 
figurations will be required to obtain the correct magnitude of downwash 
across the tailspan. 


To obtain the shape of a vortex wake behind a wing, clearly a mul- 
tiple trailing-vortex system, such as the 10-vortex systems used herein, 
is required. However, for the computation of downwash at tail locations 
8, minimum number of vortices is desirable and this number depends mainly 
on the distance between the wing wake and tail surfaces. The number 
varied from about one to three for the configurations studied herein, 
depending on whether that distance was large or small. An empirical 
relation, based on the calculations made in this study, for estimating 
the nearest integral number of vortices required per wing panel is: 


meza JA 


Epl 


where h is the height of the horizontal tail above the body axis in 
terms of bođy radii. 


Note on calculations involvi large numbers of. vortices.- In set- 
ting up the initial (trailing-edge) vortex distribution from the given 
circulation distribution I(y), one must exercise great care when n is 
large (1.6., n >20). The spacing of the vortices which approximate 
I(y) stepwise must be checked to insure a smooth curve of their divided 
differences. Otherwise, spurious effects such as the "loss" of a vortex, 
caused by excessive induced velocity, sketched below, can occur (see 
fig. 24(a)). 


"Lost" vortex 
32! 
6 5 4 C) s 58% : 


Initial vortex distribution Incorrect result at p = AL 
at u=0 due to vortex No. 2 


Sketch (x) 


a ee ee — ——— --- 


AA — ._- 


34 NACA TN 3227 


CONCLUSIONS 


An investigation has been made of a line-vortex method for computing 
vortex paths, downwash, and sidewash behind wings of wing-body combina- 
tions. Application of the method to the calculation of wake shapes, 
vortex paths, and downwash behind triangular-wing and cylindrical-body 
combinations with various prescribed data leads to the following con- 
clusions;: 


l. Wake shapes calculated by two-dimensional line-vortex theory 
agree qualitatively at least with wakes observed in & supersonic wind 
tunnel. Calculated vortex paths using this line-vortex method agree 
well with & known exact crossflow-plane solution. 


2. A missile-type wing-body-tail combination with the horizontal 
tail located in the body diametral piane lO body radii downstream of the 
wing trailing edge at 5° angle of attack was considered. It was found 
that if the plane wing is of low aspect ratio (order of 2/3), downwash 
can be computed at the horizontal tail from the wing vortices alone for 
small values of the ratio of body radius to wing semispan (order of 0.2) 
and from the body upwash alone for large values of that ratio (order 
of 0.7). For high-aspect-ratio wings (order of 4) on small bodies, down- 
wash at the tail location can be well approximated only by considering 
all the flow components - wing vortices, image vortices, and body cross- 
flow. 


3. A multiple trailing-vortex system is, of course, required to. 
determine the shape of the wake behind a wing. However, for the computa- 
tions of downwash at & tail location the number of vortices required per 
panel depends mainly on the distance between the wing wake and tail sur- 
faces. The number varies from about one to three for the configurations 
studied herein, depending on whether this distance is large or small. 

An empirical relation, based on the calculations made in this report, 
ior estimating the nearest integral number of vortices N required per 
wing panel is: 


N=1+4 Ο.16Α 

E 

— |h - μα | 

τα | 

where h is the height of the horizontal tail above the body axis, y 
the distance of the tail from the wing trailing edge (both in terms of 
body radius), α is the angle of attack in radians, A is the wing aspect' 
ratio, and r/Sm the ratio of body radius to wing semispan. 


h}. The rolling up of a vortex wake behind wings of unbanked wing- 
body combinations qualitatively resembles the wake pattern behind com- : 
parable isolated wings. The comma-shaped rolling-up pattern of the 
vortex wake behind a subsonic-leading-edge triangular wing panel, with 
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the eventual single-vortex core, is due to the nearly elliptic circula- 
tion distribution along the wing trailing edge. In contrast, a supersonic- 
leading-edge triangular wing panel generates an S-shaped vortex wake 
behind the portion of the trailing edge which has a linear circulation 
distribution, and the S-shaped wake eventually rolls up into two vortices. 


5. The spanwise distance from the body of a single vortex repre- 
senting the circulation distribution on one panel of & plane, subsonic- 
leading-edge wing and body combination is approximately constant and 
equal to about 0.76 times the exposed wing semispan for all ratios of 
body radius to wing semispan. 


6. The type of vortex wake, and resultant downwash from & wing, is 
significantly affected by the circulation distribution &long the wing 
trailing edge. It should be noted that the circulation distribution 
which must be used in setting up the initial vortex distribution &long 
the wing trailing edge is not equivalent to span loading when the wings 
are in sideslip or the body noncylindrical. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., June 1, 1954 
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NONEQUIVALENCE OF WING SPAN LOADING AND CIRCULATION DISTRIBUTION 


I(y) ALONG THE TRAILING EDGE OF A WING 


The circulation function P in fluid dynamics is defined by the 
contour integral: 


where 


ds = dxi + dyj + dzk 


Let I] be evaluated around a streamwise chord of a flat plate such 
as AB in figure 13. Then: 


c O 
regu dx, = να dx4 d uj; dx, 
O ο 


Within the limits of linear theory, u} = -uy, so that: 
ο Cu. 

r=2f Uy dx, or Esef Dam (A1) 
ο Yo o Yo 


The span loading on a wing is defined by the relation: 


wu dx, (A2) 


From the Bernoulli equation one obtains the lifting pressure coef- 
ficient for a general wing-body combination as: 


Ap -2(u-04) (v7? -vu?) + QV) (A3) 


q V yo 
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For a noncylindrical body-wing combination, two linear crossflow 
velocity potentials must be added; 9, associated with a cylindrical 
body-wing combination and 9, for an expanding body can be superposed 
as follows: 


Z 


Sketch (2) 
When the potential subscript 2 is dropped, 


uj - Uy = f(9) = -?ug » Ww, - Wy = 8(0) = "Ένα, = 0 
v2 [Yo | + [σοι Τὰ = [το [-[νφα| (Ak) 
ν΄ νη = Β|νῳ (vo, 1 + [vo D = #1 vey || va, | 


using the symmetrical properties of each potential. 


Substituting equations (Ah) into (A3), one has: 


A hu liv av 
Dee (45) 
W O O 


vhere the derivatives are evaluated on the wing (z=0). Equation (4) in 
the text is the result of substituting the expressions for potentials 
pand ® in equation (A5). It can thus be seen that the second term in 
equation (1) representing the contribution of the spanwise velocities 
(the second term of eq. (A5)), can be of the same order of magnitude as 


the first term, and not negligible as in the linearized theory for wings 
alone.t 


lThe above derivation appears in reference 21 where it is pointed 
out that for axisymmetric bodies, crossflow velocities may considerably 
exeeed axial velocities, that is, O(u) < O(v). 
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Now the span loading is obtained by substituting equation (A5) in 
(A2): 


9L qug "n Qvo 

- n dx. - η dx 

Oy af No = af Me - 
O O 


and using equation (A1): 


OL r + Ε 
zm = 2q "» = e J Vx, (Α6) 
Q O O 


Thus it is seen that the span loading along the wing trailing edge 
of a noncylindrical body-wing combination differs from the circulation 
distribution by the last integral term. In fact, it can be shown, by 
writing equation (A3) in terms of velocity potential and comparison with 
equation (7), that to any order theory span loading is not equivalent to 
circulation distribution as long as sideslip angle of flow 15 not zero 
or the body noncylindrical. For a banked cruciform configuration, there- 
fore, this nonequivalence is always present. 
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APPENDIX B 


CORRECTION OF SLENDER-BODY THEORY TO ACCOUNT 


FOR MODERATE ASPECT RATIO 


The theory developed in reference 15, for the lift of wings and 
wing-body combinations, is based on the assumption of slenderness. The 
aspect ratio of the configuration analyzed by this theory must be low 
enough so that the wings are near the axis of the Mach cone. A correc- 
tion factor is required, therefore, to apply slender-body results to | 
higher-aspect-ratio wings according to the formula: 


Lap = (κ) (TaB) opm (B1) 


The correction factor k is obtained by comparison of slender-body 
theory with the more exact linear theory (applicable to high-aspect-ratio 
wings) for triangular wings (ref. 22), that is: 


) 
k = an LI (B2) 


ui) η 


It has here been assumed that this wing aspect-ratio factor k can 
be used for the winged portion of wing-body combinations and for other 
than triangular-wing plan forms. Thus one can correct the lift of & 
higher-aspect-ratio configuration by means of equation (Bl). This assump- 
tion is justified for small values of the parameter B tan w (see fig. 26) 
when k is approximately one, and slender-body theory is itself applica- 
ble. Satisfactory results should also be obtained for large wing-to-body. 
area ratios since the wing then carries the major portion of the lift of 
the combination. Experimental data confirm the validity of these assump- 
bg wes wing-body configurations of &spect ratios on the order of 
l to 4. | í 


For triangular wings with subsonic leading edges, the lift distribu- 
tion given by linear theory has the same shape as that given by slender- 
body theory. Therefore the k factor can be defined as in equation (B2) 
by total lift ratios, and is independent of the spanwise coordinate y. 
For such wings, the method of reference 15 (applied to wings alone) yields 
results identical to the low-aspect-ratio triangular-wing results of 
Jones (ref. 23). However, it is known that Jones! results overestimate 
lift-curve slope when the parameter β tanw is not small relative to 
unity. The factor by which lift-curve slopes calculated by slender-body 
theory must be multiplied to agree with linear theory is given as: 
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pomo od 
El κ’ l-B? tan? ὦ) 
(83) 
2 
k = —— B tan w > 1 
x B tan w 


These equations have beer plotted in figure 26. 


However, for supersonic leading edges (8 tan w > 1), the lift dis- 
tributions are not functionally identical so that equation (Bl) must be 
written in terms of local chord lift. Then k becomes a function of jy; 
namely, the ratio of section lifts or circulation by linear and slender- 
body theory. Now in linear theory, the pressure coefficient on a 
supersonic-leading-edge triangular wing outboard of the Mach line is: 


Ap _ a ee | (BL) 


q 
y B^ - tan? A 


Substitution of this expression in equation (1) yields the circu- 
lation distribution (in this case equivalent to the span loading): 


Ol ο... (55) 


B^ - ten? A 


For a wing alone r - 0 in equation (3), thus giving the slender- 
body theory result: 


I(y) = 2Vga / sm? - y^ (B6) 


Hence, for the part of the wing trailing edge outboard of the Mach 
cone, the factor k for a wing-body combination is the ratio of equa- 
tion (B5) to (B6), or: 


E RENNES NEM Bm 7 Y (ΒΤ) 
/ B? tan? w-1 Sm + y 
for 
Sm “r 
r + <Y < 5m 
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Fortunately, the accuracy of the κ factor can be checked in this. 
case since an exact answer for Τ(Υ) for the wing-body combination can 
be calculated by linear theory without excessive labor. Referring to 
figure 27, at any point P(x,y) on the wing trailing edge outboard of 
the point of intersection of the Mach cone with the trailing edge, the 
potential 9, can be evaluated by integration of the distribution of 
sources at points (c,T) within the forecone of P, that is: 


p(y) = - 2 n I4 -ᾱ- dg ar (58) 


forecone 
of P 


where, in the plane of the wing (z = 0), the hyperbolic radius ry is: 


rp = y (x - 0)? - g(y-7)* (B9) 


and the angle of attack α is the sum of the body angle of attack, ap, 
and the body angle of upwash: 


e 
a = Ggli + E (B10) 


On the wing, the forecone of point P produces the sbaded area 
shown with the corresponding boundaries of integration. Thus: 


r 2 : 
V y -βίν- (+2) do 
9 (x,y)-- gp ar [^ S NIE e oe + 
3 τπτ - 
By-x T tan A A (x-o)? - βξ(γ-τ)ξ 
B-tan A 
By+x | 2 
ΠΣ | fee (a + =) do 
àT MEM ESAE 
y T ten A . X(x-c)?- g?(y-T)* 


The integrals not containing r have been evaluated in refer- 
ence 24, and in the present coordinate system give: 


(811) 


V¿ag(x-y.tan A) 


n = === c 0 (m2 
4 j p? - tan? A ( ) 


— ——— 
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The integrals containing r2/r2 yield: 


i Vopr (x-y tan A) 


fs = (B13) 
y B^y* - x* 
Therefore the total potential at P is posu s φ. or: 
Pa (X7) = - Vgag(x-y tan A) rar + ra) (B14) 
KB? -ten? A  yW Bey? - x? 


The circulation distribution T(y) along the wing trailing edge 
x = ep is then (see eq. (7)): 


"TN 
Y) ye = ody tan Asp - y) (m " ==) 
4 B^-tan? A y B°y*~ sm? tan? A 


(815) 


where 


Sm - Y 
r + 7B tan A € y € sm 


Thus equation (B15) furnishes the linear-theory answer for T(y) 
along the trailing edge outside of the Mach cone, while equation (3) 
multiplied by equation (BT) is the approximation to the linear theory 
assumed by equation (Bl) with k - k(y). This approximation is com- 
pared with linear theory in figure 28 at the spanwise position of the 
Mach line y - r + (sm-r)/ tanw for various values of r/sm as a 
function of β tan o. It is seen that the approximation assumed by 
equation (Bl) with k = k(y) is satisfactory, producing a maximum error 
of about 20 percent and only about 10 percent for the range of the 
parameter PB tan w from O to h. 


For the region of the wing inboard of the Mach cone, the pressure 
coefficient given by linear theory is (for a wing alone): 


y 
Ap ha i z p cot A - 1 — p= cot A + 1 
== sin A a Bil 
B? - tan? A p (cot A -E) x p (eot A + 4 


(B16) 
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Integration of equation (1) then results in: 


_ ΥΒ tan? w -sm 


oy 
r(y) PIDE PN. d ιο BI (Sm - y)sin7+ == -- 
E BS Eon aed * 


(Sm - y)B tan w 


(Sm + y)sin”? dl de | } (B17) 
- E 
(Sm + y)B tan w 


Therefore the factor k(y) to be used for Γ(Υ) along the wing 
trailing edge inboard of the Mach cone is the ratio of equation (B17) to 
equation (B6) or: 


= tan? ws 
k(y) = A A { om + » | (om oen? Xp o Έα, 
M (Sm? - y2)(8? tau? w - 1) (Sm -y)B tan ω 
( i yp? tan? w + Sm 
Sm + y)sin” |} 
(Sm + Υ)β tan w (B18) 
Por 
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(a) Vapor-screen photograph. 
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(b) Calculated theoretical wing vortex wake, 


Figure 1.- Vortex wake behind wing of high-aspect-ratio supersonic wing 
and body combination at angle of attack. 
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Figure 3.- Rolling up of a vortex sheet behind wing of subsonic leading-edge wing and body 
configuration 
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Figure 4.- Rolling up of a vortex sheet behind wing of supersonic lea 


Figure 5.- Rolling up of a 


Hp 


e 4 
d — 


All vortices rotating } 


A 
J ie 


ortex sheet segment containing eight equal-strength, equally spaced 
ισα. 


Lect Ni VOVN 


Ml 
MI 


NACA TN 3227 


56 
—— ——. Wing vortices 

———— —— [mage vortices 
-----------»- * Body crossflow 

ο 

o 

dD 

v 

e 

o 

ο 

i 

e 

* 

ο 

Q 


Awing tip = 5.000 
= 4045 
Wvortex 77-1350 


Avortex 


spanwise coordinate, A Sa 


(a) Aspect ratio = 2/3 


Figure 6.- Contribution of wing vortices, image vortices, and body cross- 
flow to total downwash at tail location for various aspect ratios. 
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Figure 6.- Continued. 
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Figure 7.- Vortex p&ths behind two wing-body combinations at various 


angles of attack. 
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Figure 7.- Concluded. 
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Figure 8.- Contribution of wing vortices, image vortices, and body cross- 
flow to total downwash at tail location. 
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Figure ὃ.- Continued. 
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Figure 8.- Continued. 
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Figure 8.- Continued. 
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Figure 8.- Continued. 
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Figure 9.- Vortex paths behind wings of 2/3-aspect-ratio triangular-wing 
and body combinations with various ratios of body radius to wing 
semispan r/sm- 
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Figure 10.- Variation of single-vortex strength with body radius to wing 
semispan ratio for subsonic leading-edge configuration. 
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Figure 11.- Contribution of wing vortices, image vortices, and body cross- 
flow to total downwash at tall location. 
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Figure 11.- Concluded. 
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Figure 12.- Comparison of rolling up of vortex sheet behind wing-body 
combination (A = 2/3, r/sm = 0.2) with wing alone. 
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Figure 13.- Illustration of coordinates and dimensions of general wing- 
body configuration. 
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Figure 14.- Spanwise position of single vortex approximating circulation distribution. 
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Figure 15.- Downwash along line y = z 
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Figure 17.- Calculated vortex paths behind plene wing and body combination of exemple 1. 
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Figure 18.- Calculated vortex paths behind cruciform wing end body combination of example 2. 
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Figure 19.~ Calculated vortex path behind plane wing and body combination of example 3. 
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Figure 21.- Vortex position error versus Au at stations uy = 1, 2, 3, 
and 4 for single wing vortex. 
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Figure 22.- Vortex position error.versus Au at stations 
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Figure 22.- Continued. 
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(ο) Vortex paths from p= O to ü = 4, with Ap = 1/4, 1/2, 1, and 2, viewed in crossflow plane. 


Figure 22.- Concluded. 
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Figure 23.- Comparison of stepwise approximation with Villat (exact 
crosSflow) solution for two wing-body combinations (ref. 18). 
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Figure 24.- Effect of number of vortices on downwash at tail location 
for missile-type configuration. 
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Figure 24.- Concluded. 
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Figure 25.- Effect of number of vortices on downwash at tail location 
for airplane-type configuration. 
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Figure 26.- Aspect-ratio correction factor k. 
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Figure 27.- Forecone of integration for equation (B8). 
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Figure 28.- Comparison of corrected slender-body and linear theories at Mach line on wing 
: trailing edge. 


Lect ML VOVN 


τό 


“TE 


